Ground-state properties of interacting two-component Bose gases in a hard-wall trap 
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We investigate ground-state properties of interacting two-component Bose gases in a hard-wall 
trap using both the Bethe ansatz and exact numerical diagonalization method. For equal intra- and 
inter-atomic interaction, the system is exactly solvable. Thus the exact ground state wavefunction 
and density distributions for the whole interacting regime can be obtained from the Bethe ansatz 
solutions. Since the ground state is a degenerate state with total spin S = N/2, the total density 
distribution are same for each degenerate state. The total density distribution evolves from a Gauss- 
like Bose distribution to a Fermi-like one as the repulsive interaction increases. The distribution 
of each component is N^/N of the total density distribution. This is approximately true even in 
the experimental situation. In addition the numerical results show that with the increase of inter- 
species interaction the distributions of two Tonks-Girardeau gases exhibit composite fermionization 
crossover with each component developing N peaks in the strongly interacting regime. 

PACS numbers: 67.85.-d, 67.60.Bc, 03.75.Mn 



I. INTRODUCTION 



Since two-component Bose-Einstein condensates 
(BECs) of trapped alkali atomic clouds were realized 
experimentally P, Q , low-dimensional multi-component 
Bose gases have attracted much attention from theory 
and experiment due to their connection to many areas 
of physics. Theoretical investigation mainly focuses 
on the stability, phase separation, collective excitation, 
Josephson-type oscillations and other macroscopic 
quantum many-body phenomena I3i H B 0| in the 
frame work of mean field theory. Other fundamental 
problems such as topological defects, symmetry breaking 
effects also attract growing interest [1, Q . 

Advances in experiment with ultracold atoms provide 
exciting opportunities to control and manipulate ultra- 
cold atom gases in one-dimensional (ID) waveguides by 
tightly confining the atomic cloud in two radial direc- 
tions and weakly confining it along the axial direction 
[IS [IHj [I2| • Successful realization of ID interacting quan- 
tum degenerate gases enables us to study novel many- 
body effect in various of interacting regimes, for exam- 
ple, in the strongly interacting limit, i.e. the Tonks- 
Girardeau (TG) gases 13]. Tunability of the scattering 
length cross Feshbach resonance allows experimentalists 
access to whole interaction regime from a weakly interact- 
ing limit to a strongly interacting limit. Strong correla- 
tion effect in ID quantum degenerate gases f l 4| have been 
extensively studied in recent years [l^, Il6l. llTl. [isl . [Tot . 
It is shown that ID quantum systems exhibit particular 



features which are significantly different from its three 
dimensional counterpart. 

The exact results for single component bosons with 
a repulsive (5-function interaction show that the density 
profiles evolve from Gaussian-like distribution of Bosons 
to shell-structured distribution of Fermions when inter- 
action strength increases [l^, [U, HHH, (H [H • In or- 
der to study the system in the strong interacting regime, 
non-perturbation method is highly desirable and reliable 
because the mean field theory is proven to be insuffi- 
cient. For the multi-component quantum gases, studies 
were carried out by means of various schemes such as 
mean field theory [3, HQ, extended Bose-Fermi map- 
ping in the infinitely repulsive limit f^, [13, [H, [1^ [lo] 
and exact Bethe ansatz [sil [13, [S l34| . The ID homo- 
geneous two-component Bose gas with spin-independent 
s-wave scattering (equal inter- and intra-species interact- 
ing strengthes) is exactly solvable by the Bethe ansatz 
method in a whole physical regime [sU, [s^. The fer- 
romagnetic ordering in spinor Bose gases was predicted 
some years ago 11,136 1. However, the system is no longer 
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integrable if the atomic gas no longer has the SU{2) sym- 
metry or it is trapped in an inhomogeneous potential. 
In this situation numerical methods have to be exerted 
and rich phenomena are shown for various parameters 

[2i,[13]. 

Previous study on the integrable two-component Bose 
gas mainly focuses on the energy spectrum and excitation 
properties [31i. i32l. [33l. [3^. However, important quanti- 
ties, which are related to the wave function of the system 
and are accessible experimentally, such as the density 
distribution and momentum distribution, are rarely ad- 
dressed except in the limit of infinitely repulsive inter- 
action f26|. In this paper, we are aimed to study the 
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two-component bosonic systems with SU (2) symmetry 
in the whole interacting regime by means of the Bethe 
ansatz. In the case of broken SU{2) symmetry, we resort 
to exact diagonalization method. The total density dis- 
tribution and the density distribution for each component 
can be derived from exact ground state wave function. In 
addition, numerical method will be used to evaluate the 
reduced one-body density matrix and momentum distri- 
bution of each component as well as inter- and intra- 
species density-density correlations. 

The present paper is organized as follows. Section II 
introduces the model and gives the exact solutions by 
means of the Bethe ansatz for the integrable point. In 
Section III the numerical diagonalization method is in- 
troduced and we investigate the system for accessible ex- 
perimental parameters after checking the accuracy of the 
numerical result. Section IV is devoted to the interaction 
of two TG gases. A summary is given in the last section. 



II. EXACT SOLUTION OF TWO-COMPONENT 
BOSE GAS 

We consider two-component Bose gas confined in a ID 
hard wall trap of length L, which is composed of two in- 
ternal states (pseudospin |t) and denote state 1 and 2, 
respectively) of the same kind of Bose atoms with equal 
mass mi = TO2 = "ti- The atom numbers in each compo- 
nent are A^i and N2 and N — Ni + N2 the total number. 
The many body system can be described by the second 
quantized Hamiltonian 
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a=l,2 



2ma dx dx 



+gi2 j dxii\{x)¥^{x)^2{x)^i{x), 

where ga {a = 1,2) and 512 denote the effective intra- and 
inter-species interaction which can be controlled experi- 
mentally by tuning the corresponding scattering lengthes 
ai, a2 and 012, respectively. The field operator ^J,(x) 
a{x)) creates (annihilates) an a-component boson at 
the position x. A standard rescaling procedure brings 
the Hamiltonian into a dimensionless one 



n ^ j dxY,i^i>U2 



dx"^ 



*„(x) 



+2Ui2 / dx^\{x)%{x)^2(x)^i{x) 



(1) 



with Ua — mga/fi'^ [a — 1, 2, 12). Here we have rescaled 
the energy and length in units of yi' j^vnl?' and h. The 
system with equal interaction constants Ui — U2 = 
U12 = c is integrable in both the periodic boundary 
[3l| and open boundary conditions (32l . [33| and the eigen 



problem for the original Hamiltonian is reduced to solv- 
ing the coordinate nonlinear Schrddinger equation 



with 



H^{xi,--- ,xn) =£;*(xi,--- ,xn) 



g2 



(2) 



(3) 



The Hamiltonian H commutes with the total spin oper- 
ator S, i.e., commutes with all the three components of 
the total spin operator, so that they share a common set 
of eigenstates and the system possesses a global SU(2) 
symmetry. Explicitly, the three components of total spin 
operator are defined as 

^"^l f dxY,i'lix)a'^^^,i'Ax), 

where /i, G {tii} with | (J,) corresponding to the 1th 
component (2th component) and ct" (a — x, y, z) denotes 
the Pauli matrices. 

The coordinate wave function can be determined by 
means of Bethe ansatz method and takes the following 
general form 

*(xi,--- ,xn) 

= y ^ {Xq^ — Xq^_^ ■ ■ ■ 6 {Xq2 — Xq-^ ) X 



E 



A{Q,rP) exp 



i ^ rp. kp. Xq- 



(4) 



where Q = (qi, q2, ■ ■ ■ , qn) and P = {pi,p2,- ■ ■ ,Pn) 
are one of the permutations of I,-- - ,A^, respec- 
tively, A{Q,rP) is the abbreviation of the coefficient 
Aiqi,q2,- ■ ■ ,qN;rp^Pi,rp2P2,- ■ ■ ,rp^PN) to be deter- 
mined self-consistently, and the summation (X^q) 



is done for all of them. Here r. 



± indicates that the 



particles move toward the right or the left, 0{x — y) is 
the step function and the parameters {kj} are known 
as quasi-momenta. Moreover the wavefunction Eq. 
should fulfill the the open boundary condition [37|] for 
hard wall trap 



, Xj 



0, 



^{■■- ,x,=0,---) = ^i- 
which enforces the relations 

A{Q;--- ,-Pi, ■ ■ ■ ,rp„pN) -A(Q; • • • ,pi, ■ ■ ■ ^Vp^p^) ■ 
Furthermore, the coefRcients fulfill the general relations 
A{Q;---i,j---) = Y^'A{Q;---j,i---) 



with yf given by 
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{kj ki)Pq^q^ ic 
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where Pq^qb ^he permutation operator on 

j_2f,, • • • : P) so that it interchanges qa and 
33 . Issl . l38j . For the eigenstate with total spin 
N/2 - M {0 < M < N/2), the Bethe ansatz 
equations [13, [s^ satisfied by the quasi-momentum {kj} 
and spin rapidity {Aa} are given by 



lb 
S 



exp{2ikjL) 
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n 



ki + ic kj + ki + ic 



kj iiC kj 



kj An 



ic' kj + Aq 



■A. 



id kj + Aq, + id 



N 

n 

M 
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Aa — ki — id Aa + ki — id 
Aa — ki + id Aa + ki + id 

Aa — Af^ — ic Aa + Ap — ic 
Aa — A/3 + ic Aa + Ap + ic 



with d — c/2. The energy eigenvalue is E = X^jLi kj- 
Taking the logarithm of Bethe ansatz equations, we have 



kjL 



N 

nlj — I tan 

1=1 

M 

+ ^ I tan 



— 1 kj ki 



tan 



-1 kj ~\~ ki 



1 kj Ao 



tan'i ) , (5) 



N 



A.i-v k'i 



TT Jq + I tan 



Here the quantum numbers Ij and take integer or 
half-integer values, depending on whether N — M is 
odd or even. The ground state corresponds to the case 
with M =^ [m, lal, m, [Hi. For the ground state 
Ij = {N + 1) /2—j and Ja is an empty set, and the Bethe 
ansatz equations reduce to the situation of Lieb-Liniger 
Bosegas. In this case, the coefScient ^((5;pi,P2, •■ • ,Pn) 
can be explicitly expressed as 



AiQ\Pl,P2, 

N 

i-lfUi'k 

3<l 



,Pn) 



pi ikp^ + c){ikpi + ikp^ 



with (—1)^ = ±1 denoting sign factors associated with 
even or odd permutations of P — {pi,p2, - ■ ■ ,Pn)- By 
numerically solving the sets of transcendental equations 
eq.®, the quasimomentum {kj} and thus the ground 
state wavefunction can be determined exactly. 
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FIG. 1: (color online) The ground-state density distribution of 
the first component for (a) iVi = 4, iV2 = 0, (b)iVi = 3, = 
1, (c)iVi = 2, iVa = 2, (d)iVi = 1, iVz = 3, and c = 0, 5, 15, 50. 

For the two-component Bose gas with SU(2) symmetry, 
it has been proven that the ground states are {N -1-1)- 
fold degenerate isospin 'ferromagnetic' states [sj, [H, [H, 
[36j , which are symmetrical under permutation of any two 
spins. Among the degenerate ground states, the fully 
polarized state can be represented as 



N 



N N 



J d^x * (xi, • • • xn) X 
*|(a;i)^'j(a:2)---*|(a;Ar)|0),(6) 



where ^ {xi, ■ ■ ■ xn) is given by the Bethe ansatz ground- 
state wave function (|4]). Other degenerate states can be 
generated by applying the total lowering operator S~ 
to the polarized state. For example, the total ground 
state wave function for the degenerate ground state with 
N2 spin-down particles (the state with S — N/2 and 
Sz = {Ni — N2)/2) can be expressed as 



'2' 2 



N.2 



' 2 ' 2 



where the total lowering spin operator S~ is defined as 
5- = i/ da;#|(x)*|(x) . 

In terms of the ground state wave function 
^'(xi,--- ,xjv) the total density distribution ptot(a^) = 
X]a=i 2 Pa{x) can be expressed as 



Ptot(a;) = 



iV dx2- ■ ■ dx]si\^ (x,X2, 



■dxN I* {xi,X2, ■ ■ ■ ,Xn) 



Here the ground-state density distribution of the a- 
component is given by 

Pa{x) = {^^ {S, Sz)\i'iix)^aix) {S, Sz)) . 



From the explicit form of the many body wavefunction, 
it is straightforward to get the ground-state density dis- 
tribution of the Qf-component which is found to fulfill a 
simple relation with the total density distribution (see 
Appendix) 



(8) 



We obtain unique total density profiles for all configura- 
tions with the same total atom number TV: [A''i,7V2] = 
[TV, 0] , [A^ - 1, 1] , . . . [0, TV], which is also confirmed by 
the numerical exact diagonalization method in the later 
evaluation. The conclusion ^ is valid for the integrable 
two-component boson system in the whole regime of re- 
pulsive interaction. Moreover, we would like to indi- 
cate that the conclusion ^ is valid even in the presence 
of an external confinement if the system has the total 
SU(2) symmetry, i.e., the case with Ui = U2 = U12 (soj . 
We thus recover the result in Ref. [2^ where only in- 
finitely repulsive limit was considered by a generalized 
Bose-Fermi mapping method. 

In the following calculation, L — \ will be used through 
the paper. In Fig. 1 we display the the ground-state den- 
sity distributions of the first component for A'^ = 4 and 
N2 = 0,1,2,3 for different interacting constants, where 
we find similar crossover behavior as those of a single 
component Bose gas. When the interaction is weak the 
density profiles show Gaussian-like distribution and in 
the strongly interacting regime the density profiles ex- 
hibit a shell structure with A^ peaks for each component. 
In the intermediate interacting regime the distribution 
show obvious evolution from Bose distribution to Fermi 
distribution. According to eq.®, each component takes 
the same density profiles and is normalized to the atom 
number in the component. Therefore, all configurations 
with the same total atom number share the same total 
density distribution. It is worth to note the absence of 
the demixing in the integrable system, which is contrary 
to prediction of a mean field approximation. 
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FIG. 2: (color online) The ground state energy for c = 10 
and = 4 vs utilized orbital number. Dashed line: The 
exact result of Bethe ansatz method; Scatters: Numerical 
diagonaUzation results. Units of energy: jlral? . 
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III. NUMERICAL DIAGONALIZATION 
METHOD 

Although the system is integrable for the situation of 
equal intra- and inter-atomic interactions and some exact 
results can be obtained in this case, we have to turn to 
the numerical method when the system deviates from 
the integrable point. In fact it is very difficult to adjust 
the intra- and inter-atomic interactions to be exactly the 
same in the realization of experiment. For instance the 
scattering lengths and thus the effective ID interaction 
constants are known to be in the proportion \J\ : U12 '■ 
U2 — 1.03 : 1 : 0.97 in the two components of Bose 
gas composed of internal states J^" = = —1) and 

\F — 2,mf — 1) in ®^Rb atoms In this situation 

Bethe ansatz method is not applicable and we resort to 
the numerical exact diagonalization method. 



FIG. 3: (color online) Density distribution of the ground state 
for c = 10 and A = 4. Inset: enlarged profiles in a; £ [0.3, 0.7]. 



Let us first briefly review the numerical diagonalization 
method and then investigate the ground state properties 
of the Bose-Bose mixture. The normalized eigen wave- 
function (orbital) of one particle in a hard wall takes the 

form 0i(x) — y^sin {y;x), upon which the field opera- 
tor 'i'a{x) can be expanded as 

00 

The operator (5^) creates (annihilates) one a- 
component atom in the z— th orbital. As a result the 



Hamiltonian ((T]) is discretized as 

E 

+2Ui2 ^ U,j,k,ih\J>]2^k2bii (9) 

i,j.k.l 

with — (itt)^ (i = 1, 2, 3, ...) and the dimensionless in- 
tegrals /i^-fc; = dx4>i{x)(t)j{x)(t)k{x)(t)i{x). The dimen- 
sion of the Hilbert space is now C^^^J^.J_^ x C^'^^j^^_^ if 
iVi 1-component atoms and N2 2-component atoms are 
populated on M orbitals. Then the ground state |GS) 
can be obtained after diagonalizing the Hamiltonian in 
the Hilbert space spanned by the one-particle eigenstates. 
In order to assure the precision of evaluation sufficient or- 
bitals should be considered particularly for the systems 
in strongly interacting regime. The total density distri- 
bution is given by 

2 

Ptot{x) = ^pa{x) 

with the density distribution of a-component 

p^{x) - (GSI^-Kx) *c,(a;)|GS) 

= Y,^:{x)cjy,{x){GS\h\MGS) . (10) 

In terms of the ground state wave function the reduced 
one-body density matrix for each component and the two 
body correlation of intra- and inter-species atom can be 
formulated as 

- (GS|*t(^) ^-„(x')|GS) 

= ^0*(x)0,(x')(GS|6L^.a|GS) (11) 

and 

Pal3{x,x') 

= (GS| ^l{x) ^o.{x)^\{x') ^p{x') |GS) (12) 
= ^ 0*(x)0,(a;)0:(a;')0K2;')(GS|6L&,„6l^6i^|GS). 

The momentum distribution is simply the Fourier trans- 
formation of pa{x^x')^ 

n^(k)^—[^dx[ dx'pc,{x,x')e-'<''-^'). (13) 
27r Jo Jo 

In order to test the accuracy of our numerical code, 
we compare the numerical result with that from Bethe 
ansatz method on calculating the ground state energy 
and total density profiles of two components Bose gas 
with = 4. The results are shown in Fig. 2 and Fig. 
3 for an intermediate interaction constant c = 10. The 
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FIG. 4: (color online) Density distribution of the ground state 
for TVi = 2,iV2 = 3 and U = 1.0(a), 10.0(b), 40.0(c), 80.0(d). 
PQ^(a = 1,2, tot): Bethe- Ansatz result for equal intra- and 
inter-atomic interactions {IJ\ = U2 = U12 = U); pc,{a = 
1, 2, tot): Numerical result for unequal intra- and inter-atomic 
interactions {Ui : (7i2 : U2 = 1.03 : 1 : 0.97). 



ground state energy is shown to asymptotically approach 
the BA result E/N ~ 35.22 if sufhcient orbitals are taken 
into account. For instance, we have the ground state en- 
ergy E/N = 35.60 for M — 40, the deviation of which 
is already within 1%. The density profiles calculated for 
M = 15 can match the BA result very well. In the follow- 
ing calculation, the orbital number is taken as M = 20 
(M = 15) for iV = 4 (Af = 5) and the reduced Hilbert 
space is typically composed of 10^ basis states with a 
corresponding energy cutoff (Mn)^ = 3947.84 (2220.66). 

Using the numerical method the density profiles of each 
component can be obtained even for unequal atom num- 
bers and SU{2) symmetry broken atomic interaction con- 
stants, i.e. Ui ^ U2 ^ Ui2. In Fig 4. we display both 
the total density distribution and the that of each com- 
ponent in the full interacting regime for equal and un- 
equal intra- and inter-atomic interactions in the case of 
A'^i = 2, iV2 = 3. In the situation of unique interac- 
tion constant, two components display the same density 
distribution in the full interacting regime, i.e., Pa{x) — 
^Ptot {x) ■ The density profiles show evolution from Bose 
to Fermi distribution with the increase of atomic inter- 
action. In Fig. 4d we compare the distribution of the 
system of finite strong interaction with the distribution 
of TG gas, which is obtained using the Bose- Fermi map- 
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FIG. 5: (color online) Density distribution of the ground state 
for Ni=N2 = 2 and Ui = U2 = 60. 



ping. It turns out that even if the interaction is finite the 
result from Bose-Fermi mapping can describe the system 
very well. Generally the ground state energy and the 
density profiles of two-component Bose gas do not show 
distinct difference from its single component counterpart. 
For unequal intra- and inter-atomic interacting constants 
as in the experiment {Ui : U12 : U2 ~ 1-03 : 1 : 0.97), 
the density profiles do not change drastically comparing 
with those of integrable system even in the strongly in- 
teracting regime. Particularly in the weakly interacting 
regime, the exact solution of integrable system provides 
a trustable description of the real experimental system 
because of the relatively small asymmetry of the intra- 
and inter-species interacting constants . 



IV. INTERACTION BETWEEN TWO TG 
GASES 

We have shown that how the fermionization crossover 
for the one-component Bose gas extends to a two- 
component mixture in the whole repulsive regime with 
almost equal intra- and inter-atomic interactions. The 
numerical diagonalization method can be used to inves- 
tigate two components of Bose gases even when there 
are great difference between the intra- and inter-atomic 
interactions. Now we focus on the crossover induced by 
the inter- atomic interaction constants. We start with the 
case with U12 = and Ui = U2 = 60, where each com- 
ponent is an independent TG gas, and increase the inter- 
atomic interaction U12 to see how the composite fermion- 
ization crossover happens for two initially fermionized 
components. In Fig. 5 the total density profiles are 
given for A^i = A^2 = 2 and Ui = U2 = 60. In this sit- 
uation the distributions of each component match each 
other and they are one half of the total distribution be- 
cause of equal atom numbers in two components. When 
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FIG. 6: The reduced one body density matrix for each compo- 
nent of the ground state for Ni — N2 = 2 and Ui = U2 = 60.0. 
U12 = 0.0, 30.0, 60.0 (from left to right.) Units of length: L. 
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FIG. 7: (color online) Momentum distribution of the ground 
state for iVi = iV2 = 2, Ui=U2 = 60.0. 



the inter-atomic interaction disappear the system is com- 
posed of two isolated TG gases that display peaks. 
With the increase of inter-atomic interaction the density 
profiles become flatter with more peaks appearing. In 
the strongly interacting limit the shell structure with N 
peaks display, which is the same as the density profiles 
of single component of TG gas of N atoms. For the ID 
Bose-Bose mixtures under harmonic confinement, such a 
composite-fermionization crossover has been observed as 
the interspecies cou pling st rength is varied to the limit 
of infinite repulsion [23, |30| . 

In Fig. 6 we show the reduced one body density ma- 
trix for each component, which means the probability 
that two successive measurements, one immediately fol- 
lowing the other, will find the same component parti- 
cle at the point x and a;', respectively. We notice that 
for all interacting strengthes there exists a strong en- 
hancement of the diagonal contribution pa{x,x') along 
the line x — x' . The identical momentum distributions 
for the two components are shown in Fig. 7. For all 
inter-atomic interaction strengths, Bose atoms accumu- 
late in the central regime close to zero momentum and 
the distributions decrease rapidly for large momentum. 
For strong inter-component interaction, the momentum 
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FIG. 8: Two-body correlation of intra- and inter-component 
of the ground state for Ni = N2 = 2 and Ui = U2 = 60.0. 
Top: correlation between two atoms in different components; 
bottom: correlation between two atoms in the same compo- 
nent. U12 = 0.0, 30.0, 60.0 (from left to right.) Units of length: 
L. 



in experimental, the Bose mixture shows almost same 
behaviors as the integrable system. In addition we inves- 
tigate the efibct induced by the inter-atomic interaction 
constants for two TG gases with the numerical diago- 
nalization method. It turns out that with the increase 
of interspecies interaction the system shows composite 
fermionization crossover. 
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distribution becomes broader and broader with the k = 
peak diminishing. 

It is also interesting to study the density-density cor- 
relation functions, which denote the probability that one 
measurement will find a atom at the point x and the 
other one at the point x'. In Fig. 8 we display the 
intra- and inter-species correlations between two atoms. 
At U12 = 0.0 wo have two uncorrelated TG gases and 
thus {p\{x)p2{x')) = {pi{x)) {p2{x')). With the increase 
of the inter-atomic interaction two components will try to 
avoid each other and it becomes more unlikely that one 
will find two atoms in different components at the same 
position. The intra-species correlation is always small in 
all cases because of the strong intra-atomic interactions 
in TG limit. 



V. SUMMARY 

In conclusion wc have investigated the ground state of 
two-component Bose gas with Bethe ansatz method and 
numerical diagonalization method. It turns out that the 
numerical result describe the ground state of the system 
quantitatively well. When intra- and inter-atomic inter- 
action are same (Ui = U2 = 1^12), the two-components 
Bose gas is integrable and the ground state wavefunction 
can be obtained exactly. The ground state energy and 
total density distribution are same for all configurations 
with same total atom number. With the increase of the 
interaction the total density distribution show evolution 
from a Gauss-like Bose distribution (one peak) to a shell 
structure of nointeracting spinless Fermions {N peaks). 
The distribution of each component is Na/N of the total 
density distribution. If the interaction constants deviate 
the integrable point a little, which is the real situation 



APPENDIX A 

The total ground state wave function with N2 spin- 
down particles (the state with S = N/2 and Sz = {Ni — 
N2)/2) can be generated by applying the total lowering 
operator S~ = ^ j' dxi!'^^{x)'^ ^{x) to the polarized state 
according to eq.(7) and its normalized formula shall be 
expressed as 

= j d^yi^ {xi,X2,- ■ ■ ,xn) 

x*|(a:i) • • • ¥^{xn^)^{xn^+i) ■ ■ ■ ^{xn) |0) 

with the normalized constant 
c = Ni\N2\ j dxi---dxN 

X-^ {xi,X2,- ■ ■ ,Xn)* {xi,X2,- ■ ■ ,Xn)- 

The density distribution of ath component can be ex- 
pressed as 

p^x) = (5, Sz)\ Hi^)^c.{x) 1*^ {S, Sz)) 
= 1 y d^yd^x^- (yi, • • • vn)* * (xi, ---xn) 

X (o| *T(yAr) • • • *T(yA^2+i)*i(yiV2) • • • ^i{yi)H{x) 

x*a(x)#|(a;i) • • • ^ixN2)MixN2+l) ■ ■ ■ MixN) |0) 
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For the first component the above formulation can be 
evaluated as follows 

p, (x) = (vl/^ (S, S,)\ *| {x) I M/^ {S, S,)) 

= \j c?^yd^x*(yi,---yjv)**(a;i,---a;jv) 

X (0|*T(2/Jv)---*T(2/JV2+i)*t(^) 
x^|(a;)^|(a;jv,+i) • • • ^{xn) 

^i^iiVN,) ■ ■ ■ • • • ^{xn,) |0) 

= ^iViiVilTVa! j dx2---dxN 

(x,a;2, • • • ,xn)* {x,X2, - ■ ■ ,Xn) ■ 
Similarly, for the second component we have 

P2{x) = {S,S,)\ ^{x)^ax) 1*^(5,5,)) 
= d'^yd'^x^{yi,---yN)*'b{xi,---XN) 

X (0| i'^ivN) ■ ■ ■ *T(2/JV2+i)*|(a;Ar,+i) • • • -^lixN) 
x'5^(yjvj • • • i'i(yi)i'l{x)^i{x)^{xi) ■ ■ ■ ^{xn,) |0) 

= ^N2Ni\N2\ j dx2--- dxN 

X^{X,X2,--- ,Xn)* {X,X2, - ■ ■ ,Xn)- 



Thus the total ground-state density distribution is given 

by 



Ptot{x) = Pl{x)+p2{x) 

N J dx2 ■ ■ ■ dx^^l^$ {x, X2, ■ 



,xn)* 'i' {x,X2, - ■ ■ ,Xn) 



Jdxi---dxN'b{xi,X2,--- ,Xn)* {xi,X2,--- ,Xn) 



with the density distribution of a-component 



Paix) = ^Ptot(a;). 
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